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Aims. We investigate a model for the shallow decay phases of Gamma-ray Burst (GRB) afterglows discovered by Swift/XRT in the 
first hours following a GRB event. In the context of the fireball scenario, we consider the possibility that long-lived energy injection 
from a millisecond spinning, ultramagnetic neutron star (magnetar) powers afterglow emission during this phase. 
Methods. We consider the energy evolution in a relativistic shock subject to both radiative losses and energy injection from a spinning 
down magnetar in spherical symmetry. We model the energy injection term through magnetic dipole losses and discuss an approximate 
treatment for the dynamical evolution of the blastwave. We obtain an analytic solution for the energy evolution in the shock and 
associated lightcurves. To fully illustrate the potential of our solution we calculate lightcurves for a few selected X-ray afterglows 
observed by Swift and fit them using our theoretical lightcurves. 

Results. Our solution naturally describes in a single picture the properties of the shallow decay phase and the transition to the so- 
called normal decay phase. In particular, we obtain remarkably good fits to X-ray afterglows for plausible parameters of the magnetar. 
Even though approximate, our treatment provides a step forward with respect to previously adopted approximations and provides 
additional support to the idea that a millisecond spinning (1-3 ms), ultramagnetic (B~ 10 14 - 10' 5 G) neutron star loosing spin energy 
through magnetic dipole radiation can explain the luminosity, durations and shapes of X-ray GRB afterglows. 

Key words. (Stars:) Gamma-ray burst: general — X-rays: bursts — Shock Waves — Stars: magnetars — Relativistic processes 



1. Introduction 

Before the launch of Swift in November 2004, X-ray afterglows 
of long Gamma-Ray Bursts could be pointed with X-ray tele- 
scopes not earlier than several hours after the trigger. These ob- 
servations showed in most cases a smooth power-law like decay 
F{t) oc t~ a , with typical index of a > 1. With the advent of 
Swift, X-ray fluxes could be monitored from a few minutes af- 
ter the burst trigger. These observations have revealed a complex 
behavior in the first few hours after the GRB, which nonetheless 
displays remarkably standard properties across different events. 
This behavior can be described with a double broken power law, 
with an initial very steep decay (up to few hundreds of seconds 
after the trigger) with a > 2 followed by a shallow phase, lasting 
~ 10 3 - 10 4 s, with a < 0.8 and, later on, a steeper 'normal' 
decay with a ~ 1.2-1.4 (Nousek et al. 2006, Gehrels et al. 
2009). 

The X-ray spectral slope does not change between the shal- 
low and normal decay, in marked contrast to what would be ex- 
pected in case this temporal break was caused by the passage 
of a characteristic synchrotron frequency in the X-ray band (e.g. 
Sari et al. 1998). A possible interpretation requiring no spectral 
variations in the observed energy band invokes prolonged en- 
ergy injection into the external shock that is believed to give rise 
to the GRB afterglow. Energy injection could come either from 
relativistic shells impacting the fireball at late times (e.g. Rees 
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and Meszaros 1998, Sari and Meszaros 2000) or from a long- 
lived central engine (e.g. Zhang et al. 2006, Nousek et al. 2006, 
Panaitescu et al. 2006a). 

Among different hypotheses on the nature of GRB central 
engines, two major classes can be identified. The first considers 
the formation of a black hole- debris torus system, the prompt 
emission being related to accretion of matter from the torus dur- 
ing the first ~ 10 - 100s (Narayan, Paczynski & Piran 1992, 
Woosley 1993, Meszaros, Rees & Wijers 1999). In this scenario, 
keeping the energy production active in order to power the af- 
terglow for at least ~ 10 4 s is a difficult and far from settled 
matter, (Mc Fadyen et al. 2001, Ramirez-Ruiz 2004, Cannizzo 
& Gehrels 2009, Barkov & Komissarov 2009). 

An alternative class of models invokes the formation of 
a strongly magnetic (B> 10 14 - 10 15 G) , millisecond spin- 
ning neutron star (NS, Usov 1992, Duncan & Thompson 1992, 
Blackman & Yi 1998, Kluzniak & Ruderman 1998, Wheeler 
2000). Recently, time-dependent MHD simulations have shown 
that long GRBs can originate from the interaction between a rel- 
ativistic and strongly magnetized wind produced by a newly- 
born NS and the surrounding stellar envelope. NS spin periods 
of ~ 1 ms and ultrastrong magnetic fields, i.e. B> 10 15 G, would 
be required in this case (e.g. Thompson et al. 2004, Bucciantini 
et al. 2006, 2008, 2009; see also Tchekhovskoy, McKinney & 
Narayan 2009). 

The newly formed NS is expected to loose its initial spin 
energy (> 10 52 erg) at a very high rate for the first few hours 
through magnetic dipole spin down, something that provides a 
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long-lived central engine in a very natural way. Dai & Lu (1998) 
considered this idea in relation to possible observable effects on 
the afterglow emission. Zhang & Meszaros (2001) argued that, 
in this scenario, achromatic bumps in afterglow lightcurves are 
expected for NS spin periods shorter than a few ms and magnetic 
fields stronger than several times 10 14 G. Interestingly, studies of 
the origin of NS magnetism envisage that millisecond spin pe- 
riod at birth is the key property that allows a proto-NS to amplify 
a seed magnetic field to a strength far exceeding 10 14 G, through 
efficient conversion of its initial differential rotation energy (e.g 
Duncan & Thompson 1992, Thompson & Duncan 1993). Such 
highly magnetized, fast spinning NSs are expected to loose an- 
gular momentum at a high rate in the first decades of their life 
and later become slowly rotating magnetars whose major free 
energy reservoir is in their magnetic field (Thompson & Duncan 
1995, 1996, 2001, cfr. Woods & Thompson 2006, Mereghetti 

2008) . We term these NSs as magnetars since their birth even 
though, when they spin at millisecond period, their rotational 
energy is still the main free energy reservoir. 

After the Swift discovery of early afterglow shallow phases, 
the magnetar scenario has been invoked to interpret the X-ray 
light curve of both some short and long GRBs (e.g. 05 1221 A by 
Fan and Xu 2006; 060313 by Yu and Huang 2007; GRB 050801 
by De Pasquale et al. 2007; 070110 by Troja et al. 2007). For 
GRB 060729 this scenario was shown to provide a good agree- 
ment with the shallow and normal decay phases in the optical 
and X-ray bands (Grupe et al. 2007, Xu et al. 2009). 
Finally we note that, besides the interest in understanding GRB 
physics, the very fast spin and huge magnetic field envisaged in 
the magnetar formation scenario makes these objects very in- 
teresting also for gravitational wave (GW) astronomy. Different 
possibilities for this to occur have been investigated in the litera- 
ture (Palomba 2001, Cutler 2002, Stella et al. 2005, Dall'Osso & 
Stella 2007, Dall'Osso, Shore & Stella 2009, Corsi & Meszaros 

2009) showing that, in astrophysically plausible conditions, GW 
emission might efficiently extract spin energy from the NS, in 
competition with magnetic dipole losses. The study presented in 
this paper builds on the ansatz that millisecond spinning mag- 
netars are formed in the events that give rise to long GRBs. 
We investigate the evolution of energy in a relativistic blast- 
wave subject to radiation losses due to shock deceleration in 
the ISM and energy injection from a magnetically braking NS. 
We extend previous treatments by describing the injection term 
by the standard magnetic dipole formula and deriving a predic- 
tion for the evolution of energy and luminosity that can inter- 
pret the X-ray afterglows through their shallow and normal de- 
cay phases altogether. We derive an approximate solution for the 
blastwave luminosity which we compare with X-ray GRB af- 
terglow lightcurves observed by Swift. We obtain a remarkably 
good match to these lightcurves for the range of initial spin peri- 
ods and magnetic field strengths expected for magnetars at birth. 
These results illustrate the potential of this scenario in explaining 
the early afterglow observations in a simple, unified picture. 

2. Relativistic Blast Wave with Energy Injection: 
spherically symmetric case 

We assume that a GRB event is associated to the formation of a 
millisecond spinning, ultramagnetized NS. In the context of the 
fireball scenario, the energy released in the collapse of the pro- 
genitor star produces first a fireball expanding freely at relativis- 
tic speed through the ambient medium. The prompt emission is 
produced at this early stage and is commonly ascribed to internal 
shocks in the fireball (Rees & Meszaros 1994, Paczynski & Xu 



1994, Sari & Piran 1997). A relativistic forward shock is pro- 
duced at larger distances from the explosion site (~ 10 16 cm), 
which initially propagates freely through the ambient medium. 
At a later time, call it f^, the mass swept up by the forward shock 
will be enough to begin affecting the expansion dynamics of the 
shock itself. This defines the decelaration radius » ctj, at 
which the kinetic energy of the shock starts being efficiently con- 
verted to internal energy and then radiation. This corresponds to 
the onset of the afterglow emission. We focus here only on the 
deceleration phase, describing the evolution of the total energy 
within the fireball as matter from the ISM is swept up. Our aim 
is to interpret the shallow decay phase and subsequent achro- 
matic transition to the "normal" decay phase as observed in X- 
rays, within a single physical model containing a minimal set 
of parameters. We do not address here a detailed study of the 
multiwavelength behaviour of afterglow lightcurves. In § 3.2 we 
discuss possible developments of our work in this direction, as to 
closely compare model predictions with multiwavelength obser- 
vations. The first few minutes after the GRB event are character- 
ized by a very steep power-law decay of the flux while a marked 
spectral change usually accompanies the transition to the shal- 
low decay phase (this is in contrast with the lack of spectral evo- 
lution across the shallow-to-normal transition). This initial steep 
decay is believed to arise from a different spectral component 
than the X-ray afterglow, likely the tail of the prompt emission 
(cfr. Zhang 2007 for a detailed discussion); we do not consider 
it in this work. 

In addition to deceleration in the ISM, we study the way in 
which the afterglow emission is affected by the energy injection 
caused by the spindown of the newly formed magnetar. We first 
introduce time t as that measured by a clock at rest in the NS 
(central engine) frame. In this frame the NS loses rotational en- 
ergy, likely in the form of a strongly magnetized particle wind, 
with a luminosity L«/(t) according to the usual magnetic dipole 
spindown formula 

Lsd{t) ~ (i + iKc^tf ~ (TT^F " t 2 {\ + 't/t 2 ? ' (1) 

where / is the NS moment of inertia, K = B 2 R 6 /(6Ic 3 ) with B 
the (dipole) magnetic field at the NS pole, R the NS radius and 
c the speed of light. In the second equality, the quantity L, = 
L s d(ti) represents the spindown luminosity at the initial time (t,) 
when spindown through magnetic dipole radiation sets in, and 
a = 2 Ka/j = l/f2, where t 2 represents the spindown timescale 
at time f, and w, is the initial spin frequency. E s j is the NS spin 
energy at time f„ so that L, = E S jjt 2 - The energy carried by the 
wind travels essentially at the speed of light, so that the energy 
emitted at later times by the NS can be transferred to the shock. 

To calculate the expected behavior of the lightcurve we start 
from the energy balance of the relativistic blastwave subject to 
the energy injection in eq. (1) along with radiative losses. The 
latter are described by following the prescription of Cohen, Sari 
& Piran (1998). For the time being we assume spherical sym- 
metry of all processes involved, which allows us to write the 
complete energy equation of the blast wave as 

dE E r(i) E 

— = L inj {i) - k- = (1 -B)L sd [t - — ] -k- . (2) 

at t c t 

Here k = 4e c , with e e the fraction of the total energy that is 
transferred to the electrons, r(t) is the radius of the blast wave at 
time t and all quantities are expressed in the frame of the cen- 
tral engine. Note that L,„ 7 represents the rate at which energy 
is injected in the shock at time t. This quantity is related to the 
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rate at which the central NS emitted energy - L s j - at a previous 
time, t - r(t)/c. In the central engine rest frame, an infinitesimal 
time interval dt is related to the infinitesimal displacement of the 
blast wave dr - c/3{t)dt. However, due to propagation effects, 
photons emitted at two successive radii will be received by the 
observer over a much shorter time interval dr, which defines 
what is called "the observer's time" (T). The relation between dt 
and dT is 1 



dT = (1 -p)dt 



dt 



2r 2 (?) 



(3) 



When integrated, this gives T — t - r(t)/c. Now we can trans- 
form eq. (2) to the equivalent form with respect to time T, using 
relation (3). After some manipulation one obtains 



dE 
dT 



, E I dlnf \ 
l T ldin7/ 



(4) 



In order to obtain t(T), F(f) is required (see eq. 3); this in turn 
requires a study of the hydrodynamical evolution of the blast 
wave. The solution to this problem with energy injection is far 
from trivial and beyond the scope of the present paper. Here we 
introduce an approximation in order to derive a solution to the 
problem which captures all the essential physics. Customarily 
the evolution of T in the deceleration phase is treated using self- 
similar solutions for relativistic blastwaves (Blandford & McKee 
1976, Piran 1999). In this case all quantities scale as power-laws 
(with time t). Upon writing F 2 oc r m one can solve eq. (3) and 
obtain T = t/[2(m + 1)T 2 ], from which dlnf/dlnr = 1/(1 + m) 
can be substituted in eq. (4) to obtain E(T). On the other hand, in 
the problem we are considering the shock evolution is not self- 
similar since the energy injection term introduces in the problem 
the timescale - a~ x . One can evaluate the change introduced 
by this complication using an integral expression for the total 
(internal plus kinetic) energy of the shock (Zhang & Meszaros 
2001) 



e = j/ifoipc^rV « jn(m/ 5 )rV 



(5) 



where r as ct has been assumed. By neglecting, for the sake 
of simplicity, radial variations in the density of the ambient 
medium, we can thus write F as a function of the shock energy 
and time f, namely F 2 oc E/t 3 . For self-similar solutions one ob- 
tains Eoc p~ m ; adiabatic shocks thus correspond to m = 3. The 
relation between V and E expressed by eq. (5) identifies two ex- 
tremes for the evolution of V. First, neglecting radiative losses 
correponds to the fastest possible growth rate for E which, in 
turn, corresponds to the slowest decay rate for F 2 . On the other 
hand, neglecting energy injection corresponds to the fastest pos- 
sible decay rate for E and, in turn, of F 2 . Any realistic behaviour 
of T is thus expected to lie between these two extremes. 

In the former case, which is appropriate for the early stages 
of energy injection, we can write to a very good approximation 



dE 
dt 



const 

T 2 " 



(6) 



The solution to this gives Eoc t and, thus, r oc t or dlnf/dlnr 
= 1/2. In the opposite extreme, where only radiation losses are 
present, one obtains Eoc t which implies F 2 oc r (3+k) (with, in 
general, k < 1) or dlnf/dlnr = 1/(4 + k). Note that these ex- 
tremes reproduce, as expected, the self-similar solution obtained 
by Zhang & Meszaros 2001. 



E tap50 =100 B ]4 =10 P m =l 
E„, 50 =10 B, 4 =10 P lm =l 

E = 10 B, = 3 P =1 




- 100 1000 10* 10 s 10» 10 

Observer's Time (s) 

Fig. 1. Five different theoretical (bolometric) lightcurves from the 
model presented here, drawn varying the initial energy of the afterglow 
(E; m p), the dipole magnetic field (B) and the initial spin period (P; ms ) 
of the NS. All lightcurves are obtained for the same value of k' = 0.4. 
The time at which each lightcurve begins is the deceleration time T d , 
estimated for illustration purposes simply by equating the initial energy 
Ej m p to the rest-mass energy of swept up matter in a constant density 
ISM (n =a 1 cm" 3 ). 



Hence, although the coefficient multiplying E/T in eq. (4) does 
depend on time, its value will be bracketed between the two 
extremes found above 2 , as long as £ oc r 2 f 3 holds. As these 
extremes differ by a factor 2 + k/2 at most, we can consider 
dlnf/dlnr * const as a reasonable first-order approximation, ne- 
glecting the slow and moderate change of dlnf/dlnr. This allows 
us to write the energy equation as 



dE ,E 
— =L sd (T)-k'- , 
dT T 



(7) 



where k' = ^(dlnf/dlnr)* const. Note that its value will also 
depend on the unknown density profile of the ambient medium, 
about which we do not make any assumptions. Our ignorance 
about it is completely contained in the free parameter k', as is 
our ignorance about the microphysics. In general, fixing all other 
parameters, we expect larger values of k' for a wind-like medium 
than for a constant density ISM, based on the discussion above. 
As far as our present work is concerned, the solution to eq. (7) 
can be cast in the form 



(8) 



where To > T d is any time chosen as the initial condition. The 
integral in the above expression can be expressed in terms of the 
(real valued) hypergeometric function 2F\{a, b\ c; z) 



f 



(1 +aT) 2 



dT = 



2 Fia-k',-k';2-k';z) 



,l+k 



'(k' - 1)(1 +aT) 



l-k' 



(9) 



where z = (1 + aT)~ l . Inserting this expression in the above 
eq. (8), we obtain the complete functional form of E(T) and 
can accordingly re-express the energy loss term of eq. (4) as 



the second equality holds since p ^ 1 



2 Note that, for a wind-like medium whose density declines as r~ z , 
similar conclusions would hold: adiabatic shocks correspond to m = 1, 
while dlnf/dlnr = 1 and 1/(1 + k) for the two extremes, respectively. 
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Fig. 2. Least square fit to the X-ray lightcurves of four selected GRB afterglows observed by Swift/XRT obtained through our eqs. (8) and (9). The 
blue points and red lines represent the measured (isotropic) luminosity (in the 0.3-100 keV energy range) and best fits, respectively. The starting 
time T was determined for each afterglow based on the end of the spectral transition from the previous steep-decay phase, as discussed in the text. 
Best-fit parameter values for each individual lightcurve are reported in Tab. (1). 



L(T) = k' E(T)/T, which represents the total (bolometric) lu- 
minosity of the blastwave. The resulting function thus provides 
the total bolometric luminosity as a function of (observer's) time 
T. In what follows we use it as an approximation to the X-ray 
lightcurve in order to compare it with the X-ray data. This as- 
sumes that the observed X-ray luminosity matches the bolomet- 
ric luminosity, a reasonable approximation as long as the X-ray 
flux is the dominant emission component, which is largely veri- 
fied in the early stages of the afterglow. 



3. Discussion 

In order to illustrate the salient properties of our solution, we 
plot in Fig. 1 five lightcurves calculated through eqs. (8) and (9). 
These correspond to a choice of typical values for the shock ini- 
tial energy, spin period and (dipole) magnetic field of the newly 
formed magnetar. For all curves a value k' = 0.4 for the radiative 
efficiency was assumed. 

The curves in Fig. (1) highlight some general features of the 
solution, which are not related to the nature of our approxima- 
tion. In general, the model predictions depend mainly on three 
key parameters: the NS initial spin energy, its (dipole) external 
magnetic field and the initial energy in the blast wave. The ob- 
served X-ray lightcurves of GRB afterglows display a range of 
shapes, normalizations and durations which our model can ac- 
count for in a very natural way. On average, the shallow phase is 
not flat but characterized by a negative slope, smoothly steepen- 
ing over time as energy injection decreases and radiative losses 



become more important. In fact, the NS spindown luminosity it- 
self, i.e. the energy injection term, has an effective power-law 
slope ffgd (to distinguish it from a of the observed afterglow 
lightcurve), which is expressed as Q , s d(0 = -2(f/f2)/(l + t/h)- 
This equals zero at t — and gradually steepens with time, 
reaching the value a s( j = 1 at t = t% and, eventually, — > 2 
for t —* oo. This allows a simple understanding of the lightcurve 
shapes and their dependence on the parameters. 

Given the form of the radiative loss term, ~ E/T, the condi- 
tion dE/dT > is necessary and sufficient for the shallow decay 
phase to occur, namely a part of the lightcurve where the tempo- 
ral index of the power-law decay is smaller than 1 . Therefore, the 
onset of the shallow decay phase can be defined by the require- 
ment Li n j(r) > kE(T)/T. This leads to a bound on the time at 
which the shallow phase starts, which can be expressed in terms 
of the ratio x = T/T2 as 

x E(T) 

r > k— - . (10) 

(l+x) 2 E sJ 

In general, at the initial stage of afterglow emission, the large 
value of E and small value of T make it likely that radiative 
losses dominate over the injection term. Starting from the initial 
(kinetic) energy of the explosion, £0. the afterglow luminosity 
will then decay as a power-law with index (1 + k). In this situa- 
tion, one has E(T) = Eo(To/T) k that, inserted in eq. (10) above, 
gives 



x 1+k jEotToY' 
(l+x) 2 > E,j\Tz) 



(ID 
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Putting T = 7o(<K T2) in this equation, and then defining xo = 
T0/T2, one sees that the shallow phase can start right at time To 
if jco/(1 + xq) 2 > kEo/E s j. However, since the left hand side 
in this inequality is <sc 1 by definition, the condition could be 
verified only for a very small ratio Eq/E s j. Although this might 
happen in some cases, one does not expect this to be a general 
occurrence, so that an initial power-law decay of the lightcurve, 
with index a - 1 + k > 1 is to be expected in general. After 
some time, however, condition (11) will be met so that energy 
injection will overcome radiative losses, the total energy within 
the shell will start increasing over time and the lightcurve will 
flatten accordingly. 

This shallow decay phase will clearly last as long as energy 
injection is sufficiently strong to balance radiative losses, a con- 
dition that is bound to fail somewhere after time T2. Indeed, we 
know that during the shallow phase E(T) increases over time 
oc T@, with fi < 1 and, as long as T < T2, x on the left-hand 
side of eq. (10) grows ~ T. Hence, if that inequality was sat- 
isfied at some early time, it will continue to hold at least up to 
~ T2. After time T2, on the other hand, both the left-hand side 
and right-hand side of eq. (10) start decreasing with a steepen- 
ing dependence on time. However, while the latter term has an 
asymptotic decay oc T~ k dictated by the form of radiative losses, 
the former term has a steeper asymptotic decay, oc T l , dictated 
by the form of the injection term. Therefore, the inequality ex- 
pressed in eq. (10) is bound to reverse sign at some time T > T2, 
which implies dE/dT < at that time and the afterglow lumi- 
nosity will eventually decay again oc T~( 1+k \ 

A natural anticorrelation between the duration of the shal- 
low phase and its luminosity, i.e. the afterglow luminosity, is 
also apparent from our model lightcurves, again matching one 
of the salient properties of observed X-ray afterglows (Sato et al. 
2007, Dainotti et al. 2008). In the light of the above discussion 
it is straightforward to see that this anticorrelation reflects an in- 
trinsic property of the energy injection model assumed, i.e. of 
magnetic dipole spindown. In fact, NSs with a larger spindown 
luminosity will, on average, have a shorter spindown timescale. 
This can be checked from the last step in eq. (1), in the limit 
t < t2 appropriate for the shallow decay phase, showing that the 
NS initial spindown luminosity is L, oc B 2 ^, while 1 oc B 2 oj 2 . 

This point also leads to the expectation that the most lu- 
minous among all GRBs/afterglows might even lack a shallow 
phase altogether. This expectation can also be seen directly from 
eq. (11). The left-hand side of that equation has indeed a maxi- 
mum at x = (1 + k)/(l - k) so, if the right-hand side was always 
greater than that maximum, the condition dE/dT > could 
never be met. As an example, choosing k = 0.5 gives a max- 
imum 0.325 at x = 3. Therefore, no shallow decay phase 
would occur if Eo/E s j > 0.65(72/To) 1/2 . From this example we 
see that, for a given initial spin energy, shorter values of T2 (im- 
plying stronger magnetic fields) can more easily lead to a lack 
of the shallow decay phase. On the other hand, for a given T2, 
smaller values of the initial spin energy (again implying stronger 
values of the magnetic field) favour the lack of a shallow phase. 

Finally we note that, even starting with largely different ini- 
tial luminosities, theoretical lightcurves converge to a narrow 
distribution in luminosity at late times, nicely reproducing a 
property of observed lightcurves (cfr. Nousek et al. 2006, their 
Fig. 2). Again, this can be easily understood in terms of the 
above discussion. All lightcurves evolve as r~ (1+/:) at late times, 
i.e. for t » ?2 when energy injection has long ceased, so that 
they are all parallel for a fixed value of k' . Second, the general 
anticorrelation between luminosity and duration of energy in- 
jection implies that lower luminosity plateau's will, in general, 



last longer than higher luminosity ones. Therefore the power-law 
decay starts earlier in more luminous afterglows, while less lu- 
minous ones are still (nearly) flat, which naturally causes all late 
power-law segments to look as if they had very similar normal- 
izations. Basically, this will reflect the overall energy budget of 
the blastwave, including both the initial kinetic energy and the 
energy injected in the later phase. Note that if the same total en- 
ergy is injected over a longer time - hence, at a lower luminosity 
- it can produce a late power-law decay with a slightly higher 
normalization, as clear from the second and third curves in Fig. 
1. 

3.1. Comparison with Observations and Results 

In order to further assess the goodness of our treatment, we show 
in Fig. (2) fits to four selected X-ray afterglows observed by 
Swift, based on our solution (see eq. 8). The afterglows were 
chosen for illustration purposes among GRBs with known red- 
shifts, good statistics in the XRT lightcurves, clear evidence for 
a shallow phase, sufficiently long monitoring (> 10 5 s after the 
trigger) and absence of bright flares or re-brightenings super- 
posed on their lightcurves. 

To approximate bolometric luminosity, we obtained rest 
frame 0.3-100 keV light curves (Fig. 2) from the observed 0.3-10 
keV counts rate (taken from Swift XRT lightcurve repository, see 
Evans et al. 2006, 2009), assuming an absorbed power law spec- 
tral model over the shallow phase, including also the subsequent 
normal decay if necessary. The resulting (isotropic) luminosities 
were calculated by multiplying the fluxes by 4nD 2 L , with D L the 
luminosity distance calculated by assuming a standard ACDM 
cosmology with Q. m = 0.27, Q.\ = 0.73 and Ho = 0.73 km s _1 
Mpc -1 . Finally, observer's times were corrected for the redshift 
of the source. In accordance with the idea that the initial steep 
decay phase results from the prompt event, we only fitted data 
points from the onset of the shallow decay phase, as determined 
by the time at which the spectral transition of the X-ray emission 
takes place (our To). 

We note that, even though only data points for T > To ~ 100 
s were fitted, the initial spin energy of our model refers to the 
time at which magnetic dipole braking sets in, which, in general, 
is earlier than To. Moreover, the magnetic field is assumed to re- 
main constant. We stress that our derived values of the radiative 
efficiency k' should be taken with some caution. As discussed 
in the previous section, the degeneracy between the different pa- 
rameters that determine k! can only be solved through detailed 
hydrodynamical models and joint multiwavelength fits to obser- 
vations. Here we note that best-fit values of k! correspond to val- 
ues of e e in a relatively wide, and fully acceptable, range be- 
tween 0.1 and 0.8, depending on the exact relation between / 
and T . The results of our fits show that our model is in agree- 
ment with the data in our sample. NS parameters resulting from 
our least-squares fits are reported in Tab. (1). We stress that the 
initial spin periods, ~ 1 — 3 ms, and the magnetic fields, several 
times 10 14 G up to slightly larger than 10 15 G, perfectly match 
the parameter range expected for newly born, millisecond spin- 
ning magnetars, if their seed magnetic fields were amplified to 
ultrastrong values by a strong a-Q. dynamo in the first seconds 
after formation (Duncan & Thompson 1992). In addition, we 
note that the derived magnetic field values also perfectly match 
the values of the dipole fields measured in galactic X-ray sources 
thought to host magnetars, namely the Anomalous X-ray Pulsars 
and Soft-Gamma Ray Repeaters (see Woods & Thompson 2006, 
Mereghetti 2008). Finally, the required initial spin periods imply 
spin energies in the range (3 x 10 51 - 3 x 10 52 ) erg, for a NS 
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Table 1. Best-Fit Parameters for the four selected GRBs. The time To is the time 
starting from which we fit observed lightcurves (see text for more details) and the 
energy Eo respresents the total energy, in units of 10 50 erg, within the blastwave 
at time To. The dipole magnetic field of the NS is expressed in units of 10 14 . 
Reported errors are at lcr. 



GRB 


Pi [ms] 


B14 [G] 


k 


Eo,50 [erg] 


To [s] 


050319 


1.05 ±0.03 


5.5 ± 0.4 


0.80 ± 0.25 


10.0 ±2.1 


100 


060729 


2.12 ±0.08 


3.2 ± 0.2 


0.58 ±0.11 


1.3 ±0.2 


300 


061121 


1.18 ±0.06 


12.2 ±0.7 


0.66 ± 0.21 


10.4 ± 3.0 


100 


080430 


3.7 ± 0.6 


11.7 ± 1.1 


0.33 ± 0.09 


1.2 ±0.3 


200 



moment of inertia I 1.4 x 10 45 g cm 2 (cfr. Lattimer & Prakash 
2001). Note that this is a very conservative estimate of the total 
spin energy of the NS at formation, since we are only measur- 
ing the spin energy of the NS as magnetic dipole losses set in 
as the dominant spindown mechanism. Additional mechanisms 
have been proposed, which might be initially more efficient in 
extracting the NS spin energy and angular momentum, as long as 
it spins at ^ 1 ms (cfr. Bucciantini et al. 2006, 2008, Dall'Osso, 
Stella & Shore 2009, Metzger 2010). We finally note that in at 
least one case (GRB060729), the X-ray luminosity at late times 
is somewhat lower than the model prediction. This can be caused 
by effects that have not been considered in our simple model, of 
which we breifiy mention two. First, as we discussed in § 2, 
forcing the model to fit the whole lightcurve with a fixed value 
of dlnf/dlnr can lead to an overestimate (by a factor of 2 at most) 
of the late time emission for a given shallow phase luminosity. 
This is just of the correct magnitude to explain the mismatch in 
thecaseofGRB060729. 

Secondly, fitting the X-ray lightcurves with the bolometric lu- 
minosity can lead to an overestimate, particularly at late times. 
In fact the X-ray lightcurves probably are not representative of 
the bolometric luminosity at late times, when the contribution 
of other, lower frequency bands to the total emission becomes 
non-negligible. 

3.2. Further developments 

We briefly mention here some of the points that deserve futher 
comment and that we identify as the most important steps to 
improve the model beyond the present treatment. 

As found by several authors (Panaitescu et al. 2006b, Liang 
et al. 2007), the optical lightcurves of several GRB after- 
glows show a chromatic behaviour when compared to X-ray 
lightcurves, namely they do not show any break at the shallow- 
to-normal break observed in X-rays. According to Liang et al. 
(2007), about a half of GRB afterglows with simultaneous X-ray 
and optical observations show this chromatic behaviour. This is 
customarily viewed as a major problem for most models invoked 
to explain the shallow decay phase, at least for half of the after- 
glows. As far as the model introduced here is concerned a self- 
consistent solution of T(f) is required for a detailed calculation of 
the expected synchrotron emission at different frequencies (cfr. 
Zhang & Meszaros 2004, Sari 2006). This will permit to derive 
multiband lightcurves up to late times which can be compared 
with multiwavelenght observations of GRB afterglows. However 
we stress here that a different behavior of the optical and X-ray 
lightcurves is to be expected, in general, in the synchrotron emis- 
sion from the external shock, if a spectral break frequency were 
located between the X-ray and optical bands (as already noted 
by Panaitescu et al. 2006b). This is a natural property of a flow 



in the slow cooling regime, if the cooling frequency lies above 
the optical band but below the X-ray band during (most of) the 
energy injection phase. In this case, while electrons contribut- 
ing to the X-ray emission are re-radiating instantaneously all the 
energy that is transferred to them, electrons emitting in the op- 
tical are radiating only a tiny fraction of it. Therefore, the opti- 
cal emission can evolve indipendently of the injection term for 
a while, until the emitting electrons nearly exhaust the energy 
that they accumulated. We will present a quantitative treatment 
of this issue in a future paper. 

The assumed geometry of the expanding ejecta also deserves 
some discussion. Although the subject is still open to debate, 
there exists by now considerable evidence that GRB fireballs are 
beamed, likely with opening angles of several degrees (Frail et 
al. 2001, Ghirlanda et al. 2004, Nava et al. 2007, Liang et al. 
2008, Racusin et al. 2009, Cenko et a. 2009). It is important to 
realize that the energetic requirements of the central engine as 
derived from our fits are not sensitive to the degree of beaming 
of the ejecta. This is true as long as the NS emits its spindown 
luminosity in a nearly isotropic way, as envisaged in models of 
NS magnetospheres loosing spin energy and angular momen- 
tum through the open field lines (Contopoulos, Kazanas & Fendt 
1999, Gruzinov 2005, Spitkovsky 2006). In this case, only a frac- 
tion ~ 6j of the rotational energy losses will be transferred to the 
ejecta and contribute to Lj n j in eq. (2), where 0j is the opening 
angle of the beam. On the other hand, the measured luminosity 
of the afterglow would have to be decreased by the same fac- 
tor, thus leaving the inferred (isotropic) spindown luminosity of 
the NS unaffected. The characteristic timescale f 2 of energy in- 
jection is also unaffected by beaming, implying that our derived 
values of the dipole magnetic field are also independent of beam- 
ing in this case. 

On the contrary, if the millisecond spinning NS were to emit 
its spindown luminosity in a beam of comparable size to the 
beam of the fireball, then all its power would be transferred to 
the ejecta. The required spin energy of the NS (oc w 2 ) would 
have to be reduced accordingly by a factor ~ (Pj. In this case, 

since the timescale f 2 °c B 2 coj is always independent of beam- 
ing, our derived values of the magnetic field would have to be 
increased by a factor ~ 0J 1 . 

4. Conclusions 

In the framework of prolonged energy injection models for GRB 
afterglows observed by Swift, we have considered the possibility 
that newly born magnetars - strongly magnetized and millisec- 
ond spinning NSs - are formed in the events producing (long) 
GRBs. In the first hours after formation of the NSs, the high 
spindown luminosity caused by magnetic dipole radiation losses 
represents a natural mechanism for prolonged energy injection 
in the external shock. To assess the viability of this scenario we 
considered the energy balance of a blastwave subject to injection 
of energy by a NS spinning down through magnetic dipole radi- 
ation, along with radiative losses (oc E/t). We found an approxi- 
mate expression for the (isotropic) bolometric luminosity of the 
blastwave as a function of time that is in substantial agreement 
with general properties of the shallow-decay and normal-decay 
phases of X-ray GRB afterglows observed by Swift. 

Moreover, we have shown that individual lightcurves can be 
very well fitted by using our derived expression for the bolomet- 
ric luminosity of the continuously-powered blastwave. In par- 
ticular, our best fits provide values for the initial spin period of 
the NS in the range 1-3 ms, which match well the values ex- 
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pected in magnetar formation scenarios. Best-fit values for the 
magnetic dipole field, 10 14 - 10 15 G, are also in the range ex- 
pected for such objects at formation and in agreement with the 
dipole fields estimated for Anomalous X-ray Pulsars and Soft 
Gamma-ray Repeaters, the candidate magnetars in our Galaxy. 
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